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We present an analysis of deep-UV Raman measurements of graphite, graphene and carbon nan¬ 
otubes. For excitation energies above the strong optical absorption peak at the M point in the 
Brillouin zone (« 4.7 eV), we partially suppress double-resonant scattering processes and observe 
the two-phonon density of states of carbon nanomaterials. The measured peaks are assigned to 
contributions from LO, TO, and LA phonon branches, supported by calculations of the phonon 
dispersion. Moreover, we gain access to the infrared-active Eiu mode in graphite. By lowering the 
excitation energy and thus allowing double-resonant scattering processes, we demonstrate the rise 
of the 2D mode in graphite with ultra-short phonon wave vectors. 


PACS numbers: 63.22.Rc, 78.30.-j, 61.46.Fg, 78.67.Ch 

Graphite, graphene and carbon nantoubes (CNT) have 
experienced increasing interest in fundamental research 
in the last decade. In this context, Raman spectroscopy 
has been established as a powerful experimental tech¬ 
nique, since it provides access to both the electronic and 
vibrational properties of carbon materials^’^. Due to its 
high sensitivity, it is possible to probe properties like the 
crystallographic orientation of graphene^’^, the number 
of graphene layers®, doping® and strain®’^’®, as well as 
the diameter and chiral indices (n,m) of CNT®. 

In general, the Raman spectra of graphene, graphite or 
CNT in the two-phonon region are dominated by double¬ 
resonant (DR) Raman modes®®^^®. Especially in single¬ 
layer graphene, the prominent 2D mode outperforms the 
intensity of the first-order G mode by a factor of up to 
five.^® However, these strong DR Raman modes are very 
unique only for graphitic materials. In common semi¬ 
conductors, away from optical resonances, the second- 
order spectrum typically shows the two-phonon density 
of states^^. A possible route to observe the two-phonon 
density of states also in graphene, graphite, and CNT 
is the suppression of DR Raman modes. Here, Raman 
spectroscopy with photon energies in the ultraviolet (UV) 
spectral range seems most promising, as the strong opti¬ 
cal absorption around the M point in the Brillouin zone 
is then suppressed. Thus, all so-called ’inner’ DR scat¬ 
tering processes are selectively inhibited. 

In this work, we investigate the Raman process in 
graphene, graphite, and CNT under UV excitation. For 
excitation energies well above the M-point transition 
energy of approximately 4.7 eV, we can selectively sup¬ 
press the dominant Raman processes that are commonly 
identified with ’inner’ DR Raman scattering. In these 
cases the two-phonon density of states (pDOS) is ob¬ 
served. Therefore, we gain access to phonon frequencies 
at high-symmetry points in the Brillouin zone, in par¬ 
ticular the infrared-active ifi„ mode in graphite, which 
are otherwise not accessible in either first-order or DR 
Raman scattering in the visible optical range. By lower¬ 
ing the excitation energy towards the M-point transition 



FIG. 1. Calculated phonon dispersion and the phonon density 
of states (pDOS) of graphite. Insets show the LO/TO phonon 
branches at F (Fi„ and E2g modes, which both split into LO 
and TO at finite wave vectors) and the LA, LO/TO phonon 
branches at M, respectively. Stars mark experimental results. 

energy, we can initiate inner double-resonance processes 
and therefore demonstrate the onset of the 2D mode in 
graphite. Our interpretations are supported by calcu¬ 
lations of the pDOS for graphene, graphite, and CNT. 
Furthermore, we calculate Raman spectra for graphite 
in the UV range and observe good agreement with our 
experimental data. 

We used a Horiba T64000 spectrometer equipped with 
an Argon ion laser, providing second-harmonic genera¬ 
tion of all fundamental laser lines. All spectra are cali¬ 
brated via the Raman-active vibrational mode of molec¬ 
ular nitrogen (^‘^N 2 ) from the atmosphere. We measured 
exfoliated graphite, CVD-grown graphene on Si02/Si, 
and a HiPCO-produced buckypaper CNT sample with 
a diameter range of 7 A to 13 A. In order to avoid degra¬ 
dation of our samples, we integrated Raman signals of 
graphene over a large area (80 x 80/xm^) with an inte¬ 
gration time of 80s per spectra and laser powers below 
1 mW. Integration times for graphite were in the order of 
1.5 h. Measurements were performed using a 40x deep- 
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FIG. 2. (a) Raman spectrum of graphite with an excitation energy of Sl = 5.08 eV on a logarithmic scale. Peaks are assigned 
by comparison to maxima in the pDOS (compare Fig. 1). (b) and (c) Calculated two-phonon (overtone) pDOS of 40 different 
CNT and of graphene between 2400 cm“^ and 3300 cm~^. (d) Raman spectra of a HiPCO CNT sample (black) and graphene 
(blue) with an excitation energy of 5.08 eV. 


UV objective with a numerical aperture of NA = 0.5, 
resulting in a laser spot size of ^ 1 /rm. 

Density functional theory calculations of the phonon 
spectra of graphite and graphene were performed on the 
level of the local density approximation as implemented 
in the Quantum ESPRESSO suite^®. The electrons in 
the system were modeled by projector augmented waves 
with cutoff energies of 80 Ry for the electronic wavefunc- 
tions. Reciprocal-space integration for the ground state 
was performed by grids of 21x21x1 k points for graphene 
and 21 X 21 X 6 fc points for graphite. We fully optimized 
the atomic positions and cell parameters of the consid¬ 
ered systems until the interatomic forces and the cell 
pressures were smaller than 0.001 eV/A and 0.001 GPa, 
respectively. Interactions of the slabs with residual pe¬ 
riodic images were minimized by maintaining vacuum 
layers of at least 25 A. Phonon dispersion relations were 
calculated through density functional perturbation the¬ 
ory using the same fc-point samplings as before. The 
phonon density of states was then obtained by Fourier 
interpolation of the calculated spectra onto a denser grid 
of 500 X 500 X 1 points (500 x 500 x 10 for graphite) and 
applying a Gaussian broadening of lcm“^. 

Figure 1 shows the calculated phonon dispersion of 
graphite with its density of states (pDOS). The insets 
illustrate regions with a high pDOS originating from 
the LO (longitudinal optical), TO (transverse optical), 
and the LA (longitudinal acoustic) phonon branches at 
the high symmetry points P and M, respectively. The 
Raman-active P-point vibration at ^ 1586 cm“^ (^ 2 g ir¬ 
reducible representation) is associated with the G band 
in graphene. In the case of graphite, the single-layer E 2 g 
mode gives rise to two modes with the irreducible rep¬ 
resentations E 2 g and of which the latter is infrared 
active. Its frequency is uje^^ ^ 1591 cm“^ and there¬ 
fore slightly higher than the Raman-active mode at 
~ 1586 cm“^. 

Figure 2 (a) shows a Raman spectrum of graphite at 
an excitation energy of Al = 5.08 eV. In contrast to 
previous works on UV Raman spectroscopy on graphene 


and graphite^®’we observe a distinct Raman signal in 
the two-phonon region from our samples. In fact, we 
observe hve peaks that we assign to the two-phonon den¬ 
sity of states. In detail, we observe overtone bands from 
the M point (LA/LO-derived and TO-derived phonon 
branches), the overtones of the E 2 g and modes at 
the P point and the overtone of the maximum from the 
LO-branch overbending near the P point. The infrared 
active Eiu mode in the second-order spectrum can be 
observed, since the decomposition of the direct product 
Eiu G) always contains the fully symmetric repre¬ 
sentation. Assignments are done via comparison of the 
Raman spectrum with the calculated pDOS, as indicated 
in Fig. 2(a). The intensity distribution in the pDOS 
does not fully reflect the intensity distribution in the 
measured Raman spectrum since coupling matrix ele¬ 
ments are not considered. For instance, the electron- 
phonon coupling of the LO phonon at the P point is 
large compared to LO-derived phonons with wavevec- 
tors q ^ Thus, the calculated high intensity of the 
LO-overbending pDOS is not observed experimentally. 
Instead, in the experiment, the P-point contribution is 
dominant in Fig. 2(a). We hnd the maximum of the 
LO overbending at 1626 cm“^. The frequencies of the 
Eiu and E 2 g P-point vibrations deduced from our ex¬ 
periment are 1592cm“^ and 1578cm“^, respectively, in 
good agreement with previous theoretical and experimen¬ 
tal results (Fig. 1)^®’^®. From the measured spectrum 
in Fig. 2(a), we further hnd the M-point frequency of 
the TO-derived phonon branch to be 1408 cm“^, in good 
agreement with calculations. The peak at 2688 cm“^ can¬ 
not be distinguished between the LO- and LA-derived 
branches, since they are very close in frequency at the 
M point. However, we can determine the upper bound 
of the LA-phonon branch in graphite to be 1344 cm“^. 
Thus, the phonon frequencies from the M point and the 
Eiu vibration at the P point can be investigated by op¬ 
tical spectroscopy, enabling hne adjustments for theoret¬ 
ical calculations of graphene and related systems. 

In Figure 2 (b)-(d) we show a comparison of Ra- 
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FIG. 3. (a) Equi-energy contours {Ec — = const.) from 

the electronic band structure of graphene. Highlighted is the 
contour with Ec — E^ = 5.08 eV. (b) Schematic view of a DR 
scattering process with an excitation energy of 5.08 eV. Only 
’outer’ processes can be double resonant above the transition 
energy at the M point. 


man spectra (i?L = 5.08 eV) from graphene and CNT 
with their calculated two-phonon pDOS. In contrast to 
graphite, the experimental Raman spectrum of graphene 
has a reduced complexity, i.e., less peaks are observed. 
The peaks at 3162 cm“^ and 3246 cm“^ are attributed to 
the second-order E^g T-point vibration and the second- 
order LO-branch overbending, respectively, in analogy 
to graphite. This is in good agreement with DFT cal¬ 
culations, where we find the peaks at 3154 cm“^ and 
3258cm“^. Compared to graphene and graphite, the 
UV Raman spectrum of CNTs only consists of a single 
broad, asymmetric peak around 3150cm“^. This exper¬ 
imental result is also reflected in the calculated pDOS 
for the CNT ensemble in Fig 2 (b) and can be under¬ 
stood from the fact that phonon frequencies in CNTs 
sensitively depend on the nanotube diameter and chiral 
angle^°. The broad diameter distribution of our CNT 
ensemble (TA to 13A) directly results in a broad range 
of phonon frequencies. For instance, we infer a range of 
2Aw = 80cm“^ for the second-order high-energy mode 
frequencies^^, well explaining the experimentally and 
theoretically observed broad peak centered at 3150 cm“^. 
In contrast to graphene or graphite, we do not observe 
distinct peaks in our experimental Raman spectrum that 
are related to the LO overbending or the pDOS at the 
M point. Again, also these phonon frequencies show a 
dependence on tube diameter and chiral angle^°. The 
expected range of LO-phonon maxima in our CNT sam¬ 
ple is AwLO.max = 10cm“^.^^ Thus, all contributions will 
add up to a broad shoulder on the high-frequency side of 
the main peak. The broad Raman signal towards lower 
wavenumbers is attributed to LO-, LA-, and TO-derived 
phonon bands from the M point. 

We will now turn our discussion to the DR Raman 
scattering process in the UV range. Figure 3 (a) shows 
a contour plot of equi-energy lines from the electronic 
band structure of graphene with Ec — Ey = constant {Ec- 
conduction band; valence band). We highlight the 
contour where Ec — Ey = 5.08eV, i.e., the equi-energy 
line where the absorption for an excitation wavelength 
of 244 nm occurs. In contrast to excitation in the visible 
range, UV-excitation contours form circles around the F 
point rather than closed triangles around the K point. 
Furthermore, using excitation energies above the M- 
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FIG. 4. (a) Resonant phonon wave vectors in graphite of an 
one-dimensional DR Raman process according to Eq. (1) for 
different excitation energies given next to the curves. The 
intensity for 25 eV excitation energy is multiplied by 5 • 10®. 
Resonant ’outer’ phonon wave vectors with q ~ 2.3/A cor¬ 
respond to the low-intensity Raman bands shown by dotted 
lines in (b). (b) Calculated Raman spectra according to Eq. 
(1). Dotted lines show parts of the spectra multiplied by 10^. 
All spectra are normalized and vertically offset for clarity. 


point transition energy, we can selectively suppress opti¬ 
cal transitions along the K — M high-symmetry direction 
in the double-resonance process [compare Figure 3(b)]. 
These transitions are so-called ’inner’ processes, whereas 
excitations along the K — T direction with phonons from 
K — M are ’outer’ processes. However, due to low inten¬ 
sity, the latter contribute only marginally to the Raman 
spectrum in the double-resonance process^^’^^. There¬ 
fore, the UV spectra shown in Fig. 2 are dominated by 
the (non-resonant) two-phonon pDOS signal. Tuning the 
excitation energy from the deep UV to 4.7 eV, allows op¬ 
tical transitions from the K — Mdirection and thus acti¬ 
vates dominant inner processes. 

In the following, we present and discuss calculated Ra¬ 
man spectra of graphite for excitation energies between 
25eV — 4.7eV in the spectral range from 2600cm“^ to 
3300cm“^ [see Fig. 4(b)]. DR Raman spectra were cal¬ 
culated using the equation^° 
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where we assumed the matrix elements M. to be constant. 
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FIG. 5. (a) Raman spectra of graphite for different excitation energies. The evolution of a broad and dispersive peak at 
3050 cm~^ can be clearly observed for decreasing laser energies. For better visualization, all spectra are normalized to the 
overtone of the non-resonant, infrared-active Eiu mode, (b) Comparison of experimental (black) and calculated (blue) Raman 
spectra of graphite for two excitation energies. Dashed (solid) lines correspond to the TO (LO) phonon branches. 


Exi denote the energy differences between the intermedi¬ 
ate electronic states a, b, c and the initial state i. We con¬ 
sidered the six highest phonon branches (LO, TO, LA), 
which are indexed with a. The integration is performed 
along the V — K — M high-symmetry line, including the 
scattering of both electrons and holes^^. The broadening 
factor 7 was set to 120meV for all excitation energies^^. 
The electronic bands are GW corrected, in order to fit 
the experimentally observed M point energy.^^ 

The calculated spectrum in Fig. 4 (b) for an artifi¬ 
cially high excitation energy oi El = 25 eV resembles the 
pDOS of graphite, since only off-resonant contributions 
dominate the Raman spectrum. By decreasing the exci¬ 
tation energy to values of 5 eV or less, inner DR processes 
with phonons exhibiting ultra-short wave vectors are ac¬ 
tivated and thus start to dominate the spectra [Fig. 4 (a)]. 
They mediate DR scattering processes between two elec¬ 
tronic states close to the M point. Due to the flat slope 
of the electronic bands around the M point, the double¬ 
resonance process can be fulfilled by a broad range of 
phonon wave vectors. This directly results in broad peaks 
in the Raman spectra around 3100cm“^, as can be seen 
in Fig. 4 (b). A further decrease of the excitation energy 
results in a smaller peak width and a downshift of the DR 
Raman peak, corresponding to the slope of the electronic 
bands and TO/LO phonon branches. As can be seen, the 
relative intensity of the LO-overbending Raman peak in¬ 
creases with lower excitation energies, since the resonant 
phonon wave vector approaches the maximum of the LO 
phonon branch. The maximum exhibits a wave vector 
of 0.46/A and can be resonantly accessed by excitation 
energies marginally lower than 4.7 eV (Fig. 4). 

In Figure 4(b), we also highlight the spectral range 
that is associated with ’outer’ DR processes and ex¬ 
hibits several low-intensity Raman peaks. These peaks 
are due to ’outer’ scattering processes with LA, LO and 
TO phonons and exhibit resonant phonon wave vectors 
with q « 2.3/A [Fig. 4(a)]. Thus, we can observe 
phonons with large wave vectors close to the M point 
using excitations in the deep-UV. For El = 25 eV, the 


lineshape of the low-intensity peaks directly corresponds 
to the pDOS, whereas they transform into Lorentzian 
lineshapes for lower excitation energies, since they cor¬ 
respond to outer DR processes. In general, the contri¬ 
butions from the pDOS-related Raman bands drastically 
decrease at lower excitation energies, as the enhancement 
of the DR processes outperforms the pDOS-related peaks 
by orders of magnitude. 

The experimental Raman spectra of graphite for four 
different excitation energies are presented in Fig. 5. As 
predicted by the calculations (Fig. 4), by lowering the 
excitation energy, we clearly observe the evolution of a 
broad and dispersive peak around 3050 cm“^. This peak 
is the beginning of the 2D mode in graphite, as it is 
due to a double-resonance process at the M point with 
two TO phonons exhibiting ultra-short wave vectors. As 
predicted theoretically, we further observe an intensity 
increase of the LO-overbending Raman peak for lower 
excitation energies. In Fig. 5(b), we explicitly decom¬ 
pose the calculated spectra into contributions from the 
TO and the LO phonon branch in the double-resonance 
process for two different excitation energies. It shows the 
radical change in both the experimental and calculated 
Raman spectra mainly caused by the TO phonon 
branches (blue, dashed lines), when inner processes are 
activated. This fact emphasizes the importance of inner 
contributions in DR Raman processes in graphite and 
graphene. 

In conclusion, we performed an analysis of deep-UV 
Raman spectra for graphite, graphene and CNT. At an 
excitation energy of 5.08 eV, we observe the phonon den¬ 
sity of states in the second-order spectra of all samples. 
Decreasing the excitation energy towards the transition 
energy of the M point, we can activate inner DR pro¬ 
cesses that are mediated by ultra-short phonon wave vec¬ 
tors. We observe the evolution of a broad and dispersive 
peak for excitation energies below 4.9 eV, which is the 
rise of the 2D mode in graphite. Our calculated Raman 
spectra show good agreement with the experimental data 
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and support our interpretation of the Raman process in 
the UV spectral range. Finally, our experimental findings 
can be helpful to fine-adjust calculations of the phonon 
dispersions of graphene and graphite. 
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